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We explore the results of Coffman et.al. [Phys. Rev. A, 61, 052306 (2000)] derived for general 
tripartite states in a dynamical context. We study a class of physically motivated tripartite sys- 
tems. We show that whenever entanglement sudden death occurs in one of the partitions residual 
entanglement will appear. For fourpartite systems however, the appearance of residual entangle- 
ment is not conditioned by sudden death of entanglement. We can only say that if sudden death of 
entanglement occurs in some partition there will certainly be residual entanglement. 

PACS numbers: 03.67.-a, 03.65.Yz, 03.65.Ud, 03.65.Ta 

Entanglement, a property at the heart of Quantum Mechanics, has first been brought to scientific debate the 
intriguing questions posed by Einstein, Podolsky, and Rosen in ref. [l| and since then the matter has always been 
under investigation. Recently the interest of the physical community in this counterintuitive property has raised even 
more due to its potential as a resource for information processing and quantum computation 2]. For that purpose 
having a profound knowledge of entanglement is a must (see, e.g., ref. (3] and references therein), as well as a thorough 
comprehension of entanglement distribution in cornposite systems (involving more than two degrees of freedom) . In 
this context, several years ago Coffman et. al. [J] studied the entanglement distribution in three qubit systems 
(ABC), where each one of them can be entangled with the other two. Moreover they proved the existence of what 
that they called residual entanglement, which is not detected by usual two qubits entanglement quantifiers 5, 6]. 
Their result is valid for pure states in a Hilbert space 2^2(^)2, where they proved that quantum correlation between 
A and BC will be manifest in one of three forms: i) A is entangled with B; ii) A with C; and in) the entanglement is 
distributed among ABC, the so-called residual entanglement. To this day this relation is the most general available 
in the field of quantum information. In spite of its mathematical rigor their relation has not yet been explored in 
dynamical situations. We know that entanglement distribution is very important for the implementation of quantum 
communication in general, where the relevance of entanglement distribution is crucial. This work is devoted to the 
purpose of understanding as deeply as possible relevant dynamical consequences of the relation derived in J^. 

As the work of ref. |J] was developed, an apparently disconnected effect about entanglement has been found by 
Zyczkowski et. al. Q. They very recently showed that two parties entanglement can suddenly disappear. Since then 
this dynamical characteristic of entanglement has been called sudden death of entanglement [7|-|9[ (hereafter ESD) 
and has been measured [IQi] using twin photons. A step forward in the solution to this question was given in an 
example studied by Sainz et. al. in ref. [11]. They studied a four qubit system which interacts locally and pairwise 
and showed the existence of an entanglement invariant. However in such systems entanglement sudden death is also 
present (noted first in ref. [l^). What happens to the entanglement in a unitary evolution in such a situation? Their 
result may point to the idea that the amount of quantum correlations present in a closed system should be conserved, 
however there is nothing to prevent a dynamical redistribution of the initial entanglement. In others words, the initial 
entanglement might migrate from one partition to others in a way that the initial entanglement be conserved. 

The purpose of the present work is to answer the following question: what happens with the entanglement distri- 
bution when 3-qubit systems undergo ESD? We show that residual entanglement is intimately related to ESD for a 
large class of states. 

Fourpartite systems are also investigated having the same question in mind but no solid mathematical results to 
back up our model result about the connection between ESD and appearance of genuine entanglement. In this case 
environmental effects are taken in to account. 

The context of quantum optics the kind of interaction we use and our modeling of reservoir effects has proven very 
realistic in many situations of physical interest. Since entanglement dynamics is an essential part of the implementation 
of quantum communication we believe the results presented here may be of use. 
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Tripartite-systems: Entanglement vs. Sudden Death 

Let us consider a three qubit system A, B, C where entanglement can be found in all partitions. Coffman et. al [3| 
proved that quantum correlation between A and BC will be manifest as follows 

Ca(bc)=^ab + (^ac + '^abc (1) 

where tabc stands for a tripartite residual entanglement and C^j-^ is the concurrence between partitions i and j. 
Moreover the authors noticed that tabc is invariant if one interchanges A and B. From the generality of eq.(Il]) all 
entangled physical systems which may be mapped onto a three qubits problem must obey ([T]). 

A 

I ; 13 

FIG. 1: System ABC is initially in the pure state |-i/;o) ~ \AB)\C), with AB initially entangled an factored of C. when a 
unitary interaction between A and C is turned on there will be an entanglement dynamics in the tripartite system. 

Initially for the physical situation depicted in fig. [T]we have 

Ca{bc) = Cab = C'o (2) 

and 

CaC — TABC = . (3) 

Let us consider now that the qubits A and C interact. When the interaction is "turned on" A and C will dynamically 
entangle and, according to monogomy of entanglement 13], A will be less entangled with B. However, C will "see" the 
state A as 

PA^irB{\AB){AB\) . (4) 

As during this time evolution the partitions AB and C interact one should expect that the entanglement distribution 
will be such that Cab will become smaller as a function of time and both Cac f^nd Cbc start to grow accordingly. 
This is an example where eq.(IT]) must be obeyed all along the dynamics. Therefore it is possible that besides Cab, 
Cac and Cbc there may at some point appear a tabc- So, for systems where the partition AC admits interaction 
among its constituents and shares Cq with B there are actually two very enlightening dynamical situations: i) there 
is no ESD in any of the partitions and ii) there is ESD in at least one of the partitions. In the first case (i), with 
excitation exchange between A and C, one can show that 



/o2 I /o2 

'^AB + '^BC — '-'0 



Besides this result we also have that 



must be obeyed so that 



Cb(AC) = ^AB + C%c + TABC (5) 



^B(AC) - ^0 + '^ABC ■ 

Now, since by hypothesis there is no ESD in any of the partitions and B does not interact with the partition AC ^ we 
have (^3(^(7) = Cq and tabc = 0. What happens when one of the partitions undergoes ESD (ii)? In this case, given 
the interaction between A and C the entanglement {Cac) will not disappear suddenly. Therefore ESD can only occur 
in partitions AB and BC . Let us first consider that during a time interval ESD occurs in partition AB. During this 
time window, eq.(IT]) gives 

^\(BC) = ^Ic + TABC ■ 
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However, we should remark that C "sees" ^4 as a mixed state and its capacity to entangle with A will depend on how 
much A is entangled with B and also on the type of interaction. Since initially we have C^(bc) — ^^'^ knowing 
that Cj^f^gQ-^ > Cj^fj during the whole evolution, in the interval when Cab = the residual entanglement tabc must 
be different from zero otherwise eq.(IT]) will not be satisfied. The same analysis is valid when ESD occurs in the 
partition BC, from analyzing ([S])- Last but not least we consider the case in which ESD occurs in both partitions. 
Then Cb{ac) = C'q = tabc- 

The above considerations leave no doubt that the appearence of entanglement sudden death [7|, |8[ in tripartite 
systems bears very intimate connection with higher order entanglement, i.e., residual entanglement. 

A concrete example is the tripartite system studied in ref. [l5|, consisting of two atoms, only one of which A, say, 
interacts with the cavity, the other B serves the unique purpose of allowing for an entan gled initial state with A. The 
cavity C interacts resonantly with A according to the usual Jaynes-Cummings model [16l |. where the interaction is 
given as 

Hi ^ hg{c^a^ + ca^) , (6) 

where g is a coupling constant, c (c^) is an operator that annihilates (creates) an excitation in C and = I i)(t I 
(ct+ = I t)(i I) analogously for the atoms. Consider the atomic initial state as given by 

\AB)^^=/3\n)+a\it) (7) 
and the cavity in vacuum |C)o = |0), with + jap = 1. For this initial state \ AB)^\C)o the evolved state will be 

\ABC)i"^ ^[l3cosigt)\n)+a\it)]\0)-^Psinigt)\ii)\l) . (8) 

To quantify the entanglement between A and BC of state ([5]), we will use concurrence in the form 2-\/det pA 3, where 
PA = trBc{\ABC){ABC\t) ■ For the entanglement between A and B and between A and C we use the concurrence 
which is defined in the refs. @| as 



Cp = max{0, VAi- VA2- VA3- VA4} (9) 

where the A^'s are the eigenvalues of the p {ay ® Uy) p* {ay (8) ay) organized in decreasing order, ay is one of the Pauli 
matrices and p* is the complex conjugate of p. For the state ([8]) each bipartition will have concurrence 

Cab = Col cos I 
Cac = 1/31'! sin 25^1 
Cbc = Col sin 

with Co = 2\f3a\ stands for the initial entanglement of AB and the entanglement between A and BC will be 

Ca(bc) = 2^\l3\^cos^gt{\a\^ + \l3\^sm^gt) . 
For this initial state ^ there will be no ESD in AB and AC, since 

and according to eq.® 

TABC = , 

as discussed above. Otherwise we will also have 

^A(BC} — ^AB + ^AC (10) 

showing explicitly that in this example tabc = 0- 

Next we consider an initial state which will dynamically be lead to ESD in some partition. This will happen, e.g., 
if the cavity contains one excitation initially, |C)i = For this initial state the evolved state will be 

\ABC)['^ = [l3cos{V2gt)\ n) +«cos(gi)l It)] |1) - *[asin(5t)| tt) \Q) + /3 sm{V2gt)\ U) |2)] (11) 
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FIG. 2: Graphics for the state ((TT| with P = and a = l/VTO. Left Figure (LF): Concurrences in AB and AC and 

between A and BC. The blue, red and black curves are the concurrences Cab, Cac and Ca(sc)) respectively. Right Figure 
(RF): Here we show the residual entanglement and the concurrences squared between A and B, A and C, A and BC. The 
green, blue, red and black curves are the residual entanglement tabc and the concurrences squared C\b, Cac and C^(sc), 
respectively. 



and the concurrences in AB and AC will be 

Cab — Cq max {O, I cos{gt) cos{V2gt)\ - \ sin(gt) sm{V2gt)\} 



Cac = 



\a\^\smi2gt)\ 



sm 



{2V2gt)\ 



It may be noted that ESD will be in the partition AB, as shown in Figure [5J The entanglement between A and BC 



IS 



Ca{bc) ^ 2 



y (|a|2 &in^{gt) + |/3|2 cos2(\/25i)) (\a\^ cos^{gt) + |/3|2 sin2(V2(7t)) 



and can not be written as ([TU]). It is immediate that C^c^p-, > C'^b + ^Ac ^^'^ 



^A{BC) — ^AB 
TABC > 



(12) 



(13) 



This inequality reflects the main objective of this work. Observing figure [2l we note that the residual entanglement 
exists right before of the ESD between A and B. Our interpretation of this result is as follows: the quantum 
correlations between A and B disappear for a time interval and are distributed throughout the system contributing 
to the residual entanglement [33 |. 



Residual Entanglement and Sudden Death: a conjecture 



We now consider a four partite system, A, B, C and D where initial entanglement Co is in the partition AB. We 
also consider that A interacts locally with C, and B with D as shown in figure [31 More concretely we consider two 
atoms A and B sharing an entanglement Co and the partition C consists of N oscillators initially in vacuum, same 
for D. The local interaction in the partition AC will be described by the hamiltonian 

. N N 

Hac = ^cr^ + fi'^oJkclck + h^gk{cla^ + cua^) (14) 

fe=l k=l 

where gk is a coupling constant between the atom and the k-th. oscillator of C, Ck {c\.) is the operator which annihilates 
(creates) one excitation in the fc-th oscillator of C and cr^ = | |)(t | (u^ — | t)(i I) in the atom. Similarly for BD. 

The fundamental state of the system AC, \ i)n^i does not evolve in time. However the initial state 
containing one excitation, | f) H^i \^k), evolves to the state 



l7(0>=?Wlt>|o>+x(t)U)|i) 



(15) 



FIG. 3: Initially the partition AB has Co of entanglement and is factored from CD. After an initial time a local interaction 
in the partitions AC and BD begins and there will be an entanglement dynamics in the system. 



where ^(t) and x(t) are functions to be determined which depend on N. We define the collective states 

N 

i6) = n lOfc) (16) 

k=l 

N 

ii) = {i/x{t))Y.x,mk) (17) 

k=l 

with |x(t)P = J2k=i l-^feWP and |^(i)p +Mi)P = 1 ■ When iV = 1, in the resonant limit, we have in AC and BD 
the so called Double Jaynes-Cummings [111, [ij]- The explicit forms for ^{t) and are 

m = cos(gt) (18) 
X{t) = -isinigt) . (19) 

Otherwise, when N — > oo the subsystem C is a reservoir in vacuum and A will decay exponentially as studied in 
references [1, [13, [13 ■ In this case we have 

e(t) ^ e-^*/2 (20) 
Xit) VI -e-^* (21) 

where 7 is a damping constant. 

Now let us consider the atoms prepared, as before, in \AB)ji, and the 2N oscillators in vacuum. This initial state 
dynamically evolves to 

\ABDC)['''^ ^ f3\j{t))Ac\iO)BD+a\ib)Ach{t))BD • (22) 
The concurrencies of each pair are given by 

(23) 
(24) 
(25) 
(26) 
(27) 
(28) 

From eas. (P51 - [^5)1 we may check that there will be no sudden death in any of the partitions of ABCD. The 
concurrences between A and the rest of the system is given by 



Cab 


= c,m? 


Cac 


= mMt)x{t)\ 


Cad 


= co\mx{t)\ 


Cbc 


= c^\mx{t)\ 


Cbd 


= 2\amt)x{t)\ 


CcD 


= Co\x{t)? ■ 



Caibcd) = m{t)\ VW+WbdW (29) 

which can be rewritten as 



Ca(bcd) ~ y Cab + ^ac + ^ad ' (30) 

i.e., for the atoms initially prepared in the state \AB).fp and the 2N oscillators in their vacuum state, the entanglement 
that A shares which the rest of the system is completely distributed in the partitions AB, AC, and AD. Therefore 
there will be no residual entanglement. 
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A qualitatively different situation arises if one considers the initial state 

\AB)^=(3\n)+a\U) 
for the atoms and the 2N oscillators in vacuum. This initial state evolves to the state 

\ABDC)['^^ = (3\-f{t))Ac\l{t))BD + «| ; 0)ac\ i 0)bd ■ 



(31) 



(32) 



It is well known that this initial condition, for = 1 [ll|, [l^] and — >■ oo [13, , presents ESD in some partition 
when 1/3 1 > 2\a\. We focus our attention on the entanglement that A shares with the rest of the system. The 
entanglement between A and any other subsystem and that of A with BCD are given by 



Cab = 2|/?^(t)2|max{0,|a|-|/3xW'|} 

Cad = m{t)x{t)\max{0,\a\-mt}x{m 

cac = m'mx{t)\ 

Ca(bcd) 



2mt)wmt)? + H 



(33) 
(34) 
(35) 
(36) 



It becomes apparent that the entanglement in partitions AB and AD may disappear suddenly. In the partition AB 
there will be ESD for times such that 



\x{t)? > 



In AD, ESD will occur at times such that 



\xit)W^-\x{tW > 



(37) 



(38) 



where we used the fact that \£,(t)\ = ^\ — |x(i)P • When we solve the inequality in (pS]) we find 

1 



<lx(Or< o+ \ 7 



(39) 



2 V 4 

which imposes the condition |/3| > 2\a\. So when we observe the inequalities ((37|) and ((39)) it is easy to see that when 



1 



1 



<lx(Or< 2+ V4 



(40) 



with > 2|q;|, Cab — and Cad = at the same time. This will always be the case for the initial state ([3T|) with 
|/3| > 2|a|, as shown in figurelH 





FIG. 4: Concurrences as a function oi z — |x(*)l, with /3 0.905 an a 0.429 (|/3| > 2|a|) in the state (f5^ . LF: Ca(bcd) in 
blue, Cac in red, Cad in black, and Cab in green. RF: Cl, 

BCD) in blue. Cab + C\c + C'ad in red, and Eabcd in black. 
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The entanglement between A and BCD may be rewritten as 

Cl^BCD) = Clc + CMt)\' , (41) 

which is valid during the whole evolution. The relationship above shows that the entanglement shared by A with the 
rest of the system is divided in two parts, i) one which is in AC due to the local interactions between AC and ii) one 
which is spread over the rest of the system. We may now define the positive semidefinite quantity [l8l | Eabc d which 
represents the entanglement between A and BCD which cannot be accounted for by the entanglement of A with S, 
C and D separately, i.e., Eabcd = C\(bcd) ~ \P\b + + ^\d\ ■> which for our case gives 

EABCD^Cl\m?-{ClB + Clj,\ . (42) 

Fig. S] illustrates the entanglement distribution in the case of the initial condition ([32]). Note (on the LF) that 
ESD only occurs in the partions AB and AD. When 0.584 < z < 0.812, Cad = 0. However when 0.689 < z we 
will have Cab ~ 0. In this situation when 0.689 ^ 2 ^ 0.812 we will have Cab — Cad = at the same time. 
The RF illustrates the behavior of Ca{bcd)- shows a smooth behavior when it is increasing or decreasing. It 
represents all the entanglement between A and BCD including the one coming from the unitary interaction. The 
curve in red C^g + C'^^j + C\jj initially decreases due to the entanglement decrease followed by ESD in AB and 
AD. Right after that it increases since the entanglement provided by the interaction in the partition AC becomes 
quantitatively significant. Eabcd (curve in black) presents a maximum before the other graphs. This is due to 
the fact that the entanglement between AB and between AD are decreasing and the AC entanglement is not yet 
qualitatively significant. After the ESD in AB and AD the entanglement due to the AC dynamics grows, so that the 
Eabcd curves starts to decrease. 

Interestingly enough the Eabcd entanglement will be present during the whole evolution for any value of /3 and a 
in state ([5^ . This means that for the initial state \AB)^ there will always be an entanglement between A and BCD 
which cannot be accounted for by the entanglement of A with B, C, D separately. This is not true for the initial 
state \AB)ff, where we have Eabcd = and all the entanglement content between the partitions A and BCD may 
be accounted for by two partite concurrences. When > 2|q!| in the initial state ([5T|) there will be a time interval 
At, definede by eq. (^0]) . during which Cab = and Cad = 0, as discussed above. In this situation we have 

Eabcd ^C^m\^ (43) 

which represents the entanglement distributed in the whole system that cannot be accounted for by Cab, Cac, ^^nd 
Cad- 



Discussion and Conclusion 



The interaction presented here where excitations are exchanged between the atoms and the field, simulate quantum 
circuits in Quantum Optics [19]. This interaction provides for the possibility of exchange information j20j and 
also transfer of entanglement in systems like those represented here [2l|. For example in eq. (jl4p in the limit when 
A'^ — > 00 with system C initially in the vacuum state under the well know Born-Markov approximation [225 simulates 
the vacuum fluctuations responsible for the atomic exponential decay. In spite of its broad range of applicability, this 
type of interaction does not cover all phenomena in Quantum Optics. A phase coupling between A and the A'' 
oscillators is also a useful kind of environment without excitation exchange. When N — >■ 00 this dynamics leads 
to the disappearance of coherence [2^ and may also induce ESD. In such situation one should expect that the 
entanglement distribution be very similar to the genuine entanglement in tripartite systems. This phase interaction 
between atoms and fields (when A^ — > 00) is similar to the one modeled in ref. [24] where ESD is observed when two 
entangled atoms are subjected to a classical noisy environment simulated by a stochastic classical field. This results 
in phase damping of the collective and individual atomic states. 

Recent studies show the existence of ESD in systems qubits-qutrits (2 (g) 3) 0,12^ and sudden death of nonlocality 
in three qubit systems (2 (g) 2 (g) 2) [26l - [28j have also been investigated. In the first case, as also discussed here, if 
the phase reservoir interacts either locally on globally with the qubit and qutrit, then one should expect residual 
entanglement in the subsystem of interest. The second case, where a sudden disappearance of nonlocality is observed, 
requires, however a more careful analysis since there exists entangled states which do not violate Bell inequalities 
[29143^ even when they are tripartite with residual entanglement [331]. In other words the investigation between 
entanglement and nonlocality is a very promising, open area of research. 

As for the present work we have shown that tripartite systems subjected to a local interaction will exhibit a very 
close connection between ESD and residual entanglement, based on eq. ([1]). 
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A four partite system has also been investigated and the same phenomenon is observed, Next a natural conjecture is 
in order: is ESD a general mechanics through which entanglement flows from partitions involving two qubits to larger 
ones? This is our belief based on the fact that it can be rigorously demonstrated for three qubits and several examples 
involving more qubits point in the same direction. Proving this conjecture remain an open intriguing challenge. 

The authors JGPF and MCN were partially supported by the brazilian agencies FAPEMIG (grant number CEX- 
PPM-00549-09) and CNPq (Institute Nacional de Ciencia e Tecnologia em Informagao Quantica). 
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